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Abstract. We have used sensitive archival data from the Infrared Space Observatory (ISO) to make maps of the edge-on low 
SFR galaxy, NGC 5907, in 6 different MIR bands: LW2, LW5, LW6, LW7, LW8, and LW10, covering the spectrum from 6.5 to 
15.0 yum and including several narrow bands that isolate the infrared aromatic spectral features commonly referred to as PAHs. 
Most of the MIR emission is dominated by PAHs and it is likely that emission from VSGs contribute only negligibly except 
in the broad IRAS-equivalent band. The flux ratios are typical of galaxies with low SFRs or quiesent regions within galaxies 
(e.g M 83) and a very high PAH/continuum ratio is observed. The PAH emission follows the CO distribution and also shows 
some correlation within the disk with the /1850 /jm distribution. However, the PAH emission also reaches larger galactocentric 
radii than the CO and other correlations suggest that the PAHs are also more widespread. A significant new discovery is the 
presence of PAHs in the halo of the galaxy. In the narrow bands that isolate single PAH features, the emission shows structure 
similar to high latitude features seen in other galaxies in other tracers. The features extend as far as 6.5 kpc from the plane but 
scale heights of 3.5 kpc are more typical. The /H1.3//17.7 ratio also appears to increase with distance from the major axis. To 
our knowledge, this is the first time PAHs have been seen in the halo of an external galaxy. Just as significantly, they are seen 
in a low SFR galaxy, suggesting that strong SNe and winds are not necessary for these large molecules to reach high latitudes. 

Key words, galaxies: general — galaxies: individual (NGC 5907) — galaxies: halos — galaxies: ISM 



1. Introduction 

The galactic disk-halo interface is a critical region connecting 
two very different environments over a relatively small ver- 
tical distance. The prevalence, nature, and underlying physi- 
cal drivers of features which span this region and extend into 
the halo are particularly important since they play a role in 
the energy balance as well as the chemical evolution of a 
galaxy. Moreover, knowledge of the physical conditions of 
outflowing gas in nearby galaxies can provide important con- 
straints on galaxy formation, since outflows (i.e. some form 
of 'feedback') are crucial to galaxy formation scenarios (e.g. 
Dekel & Silk [T986l Navarro & White fTff^l Mo & Mao l2TMl 
Scannapieco et al. l2"0Tj2l Kay et al. l2"0"0"2l Zhao et al. I2003al 
I2003M Sommer-Larsen et al. l2"0Tj3l Monaco |2004 Marri & 
White 2003). Nearby edge-on galaxies therefore, provide an 
important local laboratory within which such flows can be ex- 
amined in some detail. 

Evidence for outflows includes optical emission line cones 
from nuclear starbursts (e.g. Ohyama et al. 120021 Veilleux 
et al. 120031 and others), Ha disk-halo filaments (Rossa & 
Dettmar 120031 . X-ray halos (W ang, Chaves, & Irwin 120031 
Strickland et al. I2004al I2004b> . and supershells/bubbles in 



HI (e.g. Spekkens, Irwin, & Saikia 2004) and other bands 
(Lipari et al. 2004 1. Dust in halos has been observed in absorp- 
tion against the background stellar continuum (Sofue 1987; 
Ichikawa et al. HTOl Alton et al. fT^l Howk & Savage[T333| 
Thompson et al. [7004 ^ although this only probes halo dust 
fairly close to the disk. More recently, however, IR, mm, and 
sub-mm emission from dust has also been detected in halos 



(Alton et al. [19981 Neininger & Dumke [T999l Brar et al. 12003 ; 
Popescu et al. 2004, Tuffs & Popescu 2004). The physical con- 
ditions in such halos are, however, not well understood, espe- 
cially those which allow the survival of dust in hot rarefied halo 
environments. Indeed, the conditions of dust survival in halos 
may be similar to those of elliptical galaxies (Tsai & Mathews 
1995 Xilouris et al.F 
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The mid-infrared (MIR) wavelength regime has the poten- 
tial to characterise the dust properties in halos. The sensitive 
ISOCAM detector (Cesarsky et al. H996i on board the Infrared 
Space Observatory (ISO) satellite (Kessleret al. 1996) provides 
necessary spatial resolution and, with an appropriate selection 
from the vast assortment of filters (from 4 to 18 fim), it is pos- 
sible to discriminate between contributing sources of dust. By 
now, almost a decade after the launch of ISO, we are beginning 
to understand the properties of the various components in the 
MIR regime and what they are tracing in galaxies. ISOCAM 
surveys of galaxies (Rigopoulou et al. 119991 Dale et al|2000; 
Laurent et al. 2000, Roussel et al. l2001al Sturm et al. I2002i 
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and detailed studies of individual galaxies (e.g. Galliano et al. 
2003 , Vogler et al. 2005 1 have brought deep insight into the 
extragalactic MIR picture. Among the most prominent of the 
MIR features in galaxies are the carriers of the unidentified in- 
frared bands, one model for these being the PAHs (polycyclic 
aromatic hydrocarbons) 1 and small hot (typically a few hun- 
dred K, see Sect. [3} grains. 

In this study, we use the very sensitive previously unpub- 
lished observations of NGC 5907 in the ISO database. This has 
allowed us to take advantage of an extraordinary number of fil- 
ters that were used to observe both along the disk as well as into 
the halo region, where we make clear detections far from the 
plane. Combining these data, we are able to delve into the MIR 
properties of the halo. This paper is organized as follows. In 
Sect. ll.ll we review the previous observations of NGC 5907, in 
Sect. |2 we discuss the observations and data reduction, Sect. [3] 
reviews the various contributors to the MIR observing bands, 
Sect.|4]presents the results for both the disk and halo regions, 
Sect. |5] provides the discussion and the conclusions are given 
in Sect.|6] 

1.1. NGC 5907 

NGC 5907 (TableGJ, at a distance of 1 1 Mpc (1 arcsec = 53.3 
pc), is one of the largest edge-on systems in the sky. Its nucleus 
is an HII region type and no nuclear radio source has been re- 
ported (Ulvestad & Ho 2002 1 in spite of some efforts to find 
one. Non-circular molecular gas motions near the nucleus have 
been interpreted in terms of a bar (Garcia-Burillo et al. 119971 
Garcia-Burillo & Guelin fr997t : however, there is no evidence 
for a bar in the near infra-red (Jarrett et al. 2003 ). The CO dis- 
tribution is centrally peaked (Dumke et al. 19971 with a distri- 
bution suggestive of rings or spiral arms (Sofue [T994l Dumke 
et al. [T997l 

NGC 5907 has been an important target in searches for faint 
stellar halos since it has no appreciable bulge and was origi- 
nally thought to be isolated. Sackett et al. 1994 first reported 
the detection of an R-band halo to a distance of 6 kpc around 
this galaxy, following which V and I band extraplanar light was 
detected by Lequeux et al. 1996 as well as J and K band halos 
by James & Casali lT^Sl and Rudy et al. nW71 

NGC 5907 is, in fact, a member of the 396th Lyon Group of 
Galaxies (LGG 396) of which there are 4 identified members: 
NGC 5866, NGC 5879, NGC 5907, and UGC 9776 (Garcia 
I1993> . These members are at large separations (Tabled and are 
unlikely to be interacting with NGC 5907 now. However, the 
more recent detection of a nearby companion, PGC 54419, at a 
projected distance of only 36.9 kpc and velocity separation of 
A V = 45 km s 1 (Shang et al. [1^981 indicate that NGC 5907 
is not as isolated as previously thought. Moreover, their addi- 
tional discovery of a large optical ring around the galaxy indi- 
cates that some interaction has indeed occurred. A pronounced 
HI warp (Sancisi 1976 1 and some evidence for an optical warp 
(Morrison et al. 1199 9) are consistent with this. The ring and 

1 The carriers of the MIR bands are not certain, but in this paper, 
we adopt the PAH nomenclature for discussion and comparison with 
other authors. 



nearby companion are faint in comparison to NGC 5907; the 
luminosity of the ring is < 1.2% of NGC 5907 and the mass 
of the second companion is 0.5% of the mass of NGC 5907 
(Shang et al. 1998 1. Dynamical modelling of the ring indicates 
that it was formed at least 0.8 Gyr ago via the destruction of a 
satellite of mass 2 x 10 8 M (Johnston et al. 2001 who adopted 
a slightly larger distance to the galaxy). Follow-up observa- 
tions showed that the ring is highly asymmetric and suggests 
that NGC 5907 does not after all have a faint extended stel- 
lar halo as originally envisioned (Zheng et al. 1999). A near- 
infrared (3.5-5 fim) search also failed to show evidence for 
a halo associated with NGC 5907 (Yost et al.[T333) and HST 
observations find fewer bright giants in the halo region than 
would be expected from a halo with standard dwarf-to-giant 
ratios (Zepf et al. 2000). Recent 2MASS K s band observations 
of NGC 5907 (Bizyaev & Mitronova l2"0"0"2T ) provide a stellar 
scale height of zo = 0.49 kpc with an isothermal assumption 
[i.e. p(z) oc sec 2 (z/zo)]. 

As for a gaseous halo, no high latitude ionized gas has yet 
been detected (Rand 1006) nor did Howk & Savage [1999] de- 
tect any extraplanar dust in optical absorption. Dumke et al. 
nWl find a CO FWHM of 400 pc and an HI FWHM of 1.5 
kpc, although the warp makes this latter value uncertain. The 
most extended component appears to be the radio continuum 
(Dumke et al. 2000), for which an exponential thick disk com- 
ponent is found with a scale height of 1.5 kpc. An 850 fj.m map 
has recently been published by Alton et al. (2004) who find an 
extraplanar exponential scale height of 0. 1 1 kpc. 

ISO 60, 100, and 180 //m fluxes as well as a 12 /urn im- 
age have been published by Bendo et al. (|2002 ), and Bendo et 
al. l|2003l also present temperature fits to the data. More re- 
cently, Alton et al. have also fit a two-component dust model 
to NGC 5907 which includes newer 850 fim data from ob- 
servations using the Submillimetre Common Use Bolometric 
Array (SCUBA) on the James Clerk Maxwell Telescope. These 
temperatures are presented in Tabled An earlier 2-component 
dust temperature model of 18 K and 54 K had been reported 
by Dumke et al. {1997) who observed the 1.2 mm dust com- 
ponent. Note that two-temperature models are simply a rough 
approximation to the range of temperatures that likely exist in 
this galaxy. Other parameters for NGC 5907 are also listed in 
Tabled including the star formation rate (SFR) which is quite 
low (cf. SFR = 3.8 M Q yr l for the quiescent galaxy, NGC 891, 
Popescu et al. 2002). In Table ^ we also provide comparative 
data on the nearby galaxy, M 83, whose spectral energy distri- 
bution (SED) has recently been studied in some detail (Vogler 
et al. 2005 ) and with which we compare our NGC 5907 results. 

2. Observations & Data Reduction 

2.1. Observations 

Data on NGC 5907 were available in filters, denoted LWn, 
where LW indicates the long wavelength array of the ISOCAM 
detector and n specifies the wavelength band. Details of the 
data acquisition and resulting maps are provided in Table [3] 
As we prefer to designate the bands according to their wave- 
length, we use the following nomenclature: LW5, LW6, LW7, 
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Table 1. Basic Galaxy Parameters 



Parameter 


NGC 5907 


M 83 


RA (J2000) a (h m s) 


15 15 53.69 




DEC (J2000) a (° ' ") 


56 19 43.9 




V/,„/ a (km s~ 1 ) 

* net \ / 


667 




Distance (Mpc) 


ll b 


4.5 C 


Morph. Type a 


SA(s)c: sp ffll: 


SAB(s)c; HII 


Incl. (°) 


86.5 d 


27 c 


Major x Minor Axis Diameters a (' x ') 


12.77 x 1.40 


12.9 x 11.5 


Major Axis Diameter (kpc) 


41 


17 


/kn / /"i nn^" 

J out J LUU 


0.24 


0.54 


IK \ U/ 


7.1 x 10' 


1.4 x 10 10 


T rf „„ 6 1-Comp M" 1 ) 


26 K 


36 


T rf „„ g 1-Comp (,t~ 2 ) 


22 K 


29 


T dus , g 2-Comp (/T 2 ) 


13, 26 K 




W 2-Comp (A- 15 ) 


13, 28 K 




M dml (Mq) 


6.0 x 10 7 \ 3.2 x 10 7i 


1.8 x 10 6 ' 


Z rf ,,„ J (Mo kpc" 2 ) 


7.9 xlO 5 


2.6 x 10 5 


SFR (M Q yr-') 


2.2\ 1.2 k 


2.4i 


SFR/A 1 (Mo yr 1 kpc" 2 ) 


1.7 x 10~ 3 , 9.0 x 10~ 4 


1.1 x 10~ 2 



a From the NASA Extragalactic Database (NED). This center is the 2MASS center which agrees with the dynamical center (Dumke et al. 

[19971 to within A RA = 0.8" and A DEC = 0.3" . 
b Sasaki [T987l 

c Thim et al. l2003l 

d Garcia-Burillo et al. [T997l 

e From the major and minor axis diameters and a thin disk assumption. 

f Flux ratio or FIR luminosity, using the 60 and 100 yum IRAS fluxes from Sanders et al. 2003 L IR allows for extrapolations below 40 and 

above 120 yum and adjusts to the distance given in this table. 
g Dust temperature, assuming the emissivity laws shown, for one or two component fits (Bendo et al. 2003 1. For M 83, a 2-component model 

does not improve the fit. The second line for the 2-component model gives the fit from Alton et al. 1 2004 1. 
h From the 2-component dust model of Misiriotis et al. 1200 II 

' From Bendo et al. 1 2003 1 for fluxes within a 135" aperture, assuming a AT 1 emissivity and adjusting to the distances in this table. 
1 Dust surface density, using the dust masses in this table (from the adjusted Bendo et al. data) and their aperture sizes. 
k From L ;R of this table and the calibration of Kennicutt l 19981 

1 SFR per unit area, from values in the immediately preceding row and the optical major axes given in this table. 
Table 2. Galaxy Group Membership 3 



Galaxy Name 


Velocity 


RA 


DEC 


Separation 15 


Separation 11 




(km s" 1 ) 


(h m s) 


C ' ") 


(Arcmin) 


(Radii) 


NGC 5907 


667 


15 15 53.69 


56 19 43.9 






NGC 5866 


672 


15 06 29.56 


55 45 47.9 


85.8 


13.4 


NGC 5879 


772 


15 09 46.77 


57 00 00.7 


64.5 


10.1 


UGC 9776 


833 


15 13 07.44 


56 58 07.3 


44.7 


7.0 


PGC 54419 


712 


15 14 48.04 


56 27 15.4 


11.8 


1.8 



a Members of the Galaxy Group LGG 396 (Garcia 1993 1 as well as PGC companion found by Shang et al. 1998 All values are taken from 
NED. 

b 'Separation' indicates the distance of the group member from the center of NGC 5907 in units of arcmin and in units of the radius of 
NGC 5907 (see Table[TJ. 



LW8, LW2 and LW10 are referred to as 6.8N, 7.7N, 9.6N, and 
11. 3N, 6.7W, and 12W, respectively, where the number repre- 
sents the central wavelength of the filter and N/W refers to a 
'narrow'/' wide' band. 

The first 4 data sets were taken on 23 Feb 1996 in stan- 
dard raster mode (CAM01) in the narrow bands, 6.8N, 7.7N, 
9.6N, and 11. 3N. Each observation had 3 different on-source 



pointings for a total field which provided full coverage of the 
south-east part of the galaxy as well as sufficient coverage far 
from the emission for sky subtraction. For each pointing, a se- 
ries of frames were taken over a period of time in a 'temporal 
block' . The approximate number and duration of the frames are 
given in Table [3] 
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Table 3. Observing & Map Parameters 



Parameter 


6.8N (LW5) 


7.7N (LW6) 


9.6N (LW7) 


11. 3N (LW8) 


12W (LW10) 


6.7W (LW2) 


Central Wavelength 11 (//m) 


6.8 


7.7 


9.6 


11.3 


12.0 


6.7 


Frequency Range 11 (//m) 


6.5 - 7.0 


7.0-8.5 


8.5 - 10.7 


10.7 - 12.0 


8.0 - 15.0 


5.0 - 8.5 


Observing Mode ab 


CAM01 


CAM01 


CAM01 


CAM01 


CAM01 


CAM03 


TDT No. b 


09801401 


09801401 


09801401 


09801401 


26902911 


15601301 














15601202 














15601103 














15601004 


Pixel Field of View (arcsec) 


6.0 


6.0 


6.0 


6.0 


1.0 


6.0 


PSF (FWHM) (arcsec) 


7.2 


7.2 


7.2 


7.2 


4.2 


7.2 


Date of Observations 


23 Feb 96 


23 Feb 96 


23 Feb 96 


23 Feb 96 


12 Aug 96 


21 Apr 96 


No. of On-Source Pointings 6 


3 


3 


3 


3 


12 


4 


Mean No. Frames per Pointing 


16 


16 


17 


17 


15 


248 


Integration Time per Frame (s) 


10.0804 


10.0804 


10.0804 


10.0804 


1.12 


10.0804 


Total On-Source Time' 1 (min) 


8.1 


8.1 


8.6 


8.6 


3.4 


167 


Sky Coverage 6 (arcmin 2 ) 


8.2 x 3.2 


8.2 x 3.2 


8.2 x 3.2 


8.2 x 3.2 


2.8 x 1.6 


8.4 


Calibration Error" (%) 


7.5 


6.4 


6.7 


5.4 


3.9 


3.3 


Sky Level f (mjy pixel -1 ) 


1.39 


3.28 


7.19 


10.27 


0.27 


1.99 


l<x Sky Error g (mjy pixel" 1 ) 


0.24 


0.15 


0.28 


0.30 


0.015 




cr Rms h (mjy pixel" 1 ) 


0.44 


0.17 


0.15 


0.24 


0.062 




Transient Rms 1 (mjy pixel" 1 ) 


0.35 


0.24 


0.16 


0.23 


0.044 




cr 1 (mjy pixel" 1 ) 


0.61 


0.33 


0.36 


0.45 


0.077 


0.025 k 


cr 1 (mjy arcsec" 2 ) 


0.017 


0.0092 


0.010 


0.013 


0.077 


0.00069 


Dynamic Range" 1 


20 


72 


28 


41 


6.9 


702 



a Blommaert et al. 120031 

b CAM01 is the general observing mode for photometric imaging. CAM03 is a beam-switching mode for photometric imaging. TDT no. is 

a unique number that identifies the observation. 
c This is the number of different pointing positions on-source; for some pointings, the fields overlap. 
d Total integration time per on-source pointing times the number of pointings. Time on sky not included. 
e Total area of sky covered after registration of the different on-source pointings. 
f Median value of the histogram of all sky pixel values which were subtracted from the image. 
g ltr width of the histogram of sky pixel values which were used to determine the sky background median. 

h Median a Rms over the field, where the <x Rm , is the per pixel error from the dispersion in the signal carried through by the software throughout 

the data reductions. 
1 Median rms over the field, where the rms is the error due to transients. 

1 Quadratic sum of the random errors from the immediately preceding 3 rows, unless otherwise indicated. This error represents a typical 
random error for any pixel in the map, but there are positions at which the error is larger or smaller. Note that the calibration error is not 
included. 

k Median rms error over the 4 fields which constitute the mosaic, including all applicable random errors. For this beam-switching observation, 

the above errors were factored into this error in the software and were not calculated separately. 
1 Random error from previous row per 1 arcsec square pixel. 
m Maximum of map divided by the median Random error. 



The 5th data set was taken on 12 Aug 1996 (12W). These 
data were taken in microscan raster mode (CAM01) with 12 
observing positions in a 6 x 2 pattern (Bendo et al. 2002 1 with a 
resulting 1" pixel spacing. In this case, 12 on-source pointings 
were made for full coverage of the central part of the galaxy 
with sufficient sky coverage for subtraction. 

The final data set, 6.7W, was obtained on 21 Apr 1996 in 
beam switching mode (CAM03). There were 4 positions on- 
source: one containing the galaxy center but offset to the south- 
west, another containing the galaxy center offset to the north- 
east, one at the end of the major axis to the north and one at 
the end of the major axis at the south. Thus, the final combined 
image does not provide full coverage of the whole galaxy but 
rather samples the galaxy at four locations. For each of the 4 



pointings, the off-source sky positions were placed at a variety 
of position angles with respect to that pointing. 

2.2. Reductions 

All data were reduced using the CAM Interactive Reduction 
Package (CIR, Chanial 2003 1. First, the dark current was sub- 
tracted following the method of Biviano et al. (1998 1 which 
includes a dark correction, a second-order dark correction de- 
pending on detector temperature and time of observation, and 
a short-drift correction. Next, high glitches due to cosmic ray 
(CR) hits were removed automatically via a multi-resolution 
filtering technique (Starck et al. 1999J at a 6<x level. The ef- 
fects of memory on a pixel (transient effects) were then cor- 
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Fig. 6. As in Fig. but for 6.7W with contours at 0.0008 
{2cr map ), 0.0012, 0.002, 0.005, 0.010, 0.030, 0.10, and 0.35 mJy 
arcsec -2 . The irregular region in (b) denotes the area common 
to all maps. 



rected (de-glitched) using the Fouks-Schubert method (Coulais 
& Abergel 2000). Each data set was then examined carefully, 
frame by frame, and further bad pixels were removed manually 
as required. Most of these occurred in regions immediately ad- 
jacent to pixels which had been automatically de-glitched. 

For the raster mode observations (CAM01: 6.8N, 7.7N, 
9.6N, 11. 3N, and 12W), the edges were then trimmed 
(blanked), the frames were corrected for jitter, and averaged. 
The results at each pointing were then corrected for the fiat 
field according to the on-line library of calibration flat fields 
(see Roussel et al. 12001 at . The various pointings on-source 
were then projected onto a common grid in sky coordinates and 
corrected for field lens distortion. At this point, the data were 
calibrated to mJy/pixel units (Blommaert et al. 2003 1. Finally, 
the sky background/foreground was subtracted by first isolat- 
ing a region of sky-only in the field and subtracting the median 
of that region from the field. 

For the beam switching observations (CAM03: 6.7W) we 
carried out an additional stage of manual deglitching to account 
for the fact that in the off-source positions, the CCD contained 
a fading memory of galaxy itself from when the array pointed 
at the galaxy. Thus, the region of the array which had previ- 
ously contained the galaxy was blanked in approximately 10 
sky frames immediately after moving from the galaxy to the 
sky. This blanking covered about 20% of the array and left ap- 
proximately 20 frames where the image had faded and blank- 
ing was no longer required. This additional deglitching was re- 
quired only for the center-most two fields, and not for the two 
fields at the ends of the galaxy's major axis, given the weak- 
ness of the signal there. Beam switching allowed for sufficient 
sky coverage off source that the sky positions provided both the 
flat field and its distortions as well as the sky brightness values. 
Thus, after the average on-source image was created, it was cor- 
rected using an image that was formed from an average of the 
various off-source positions. The result was then calibrated to 
mJy/pixel units. This mode of operation provides the best sky 
subtraction because it does this pixel by pixel rather than using 
a single value (the median of the sky values) for all pixels, as 
above. Finally the registration and combination of the 4 point- 
ings, as well as an interpolation over one bad column on the 
array was carried out outside of the CIR package. In the region 
of the galaxy in which there was overlap (see center region of 
Fig-Et), the flux levels between the two overlapping fields was 
~ 4.5%. 

All images were interpolated/rotated into RA/DEC coordi- 
nates with 1 arcsec square pixels for consistency. The resulting 
maps are shown in Fig.^ through^. The lowest contour, in 
each case, is set to 2cr mnp , where (T map is the rms variation cal- 
culated from a blank sky region on the emission map itself. 
The field of view is accurately shown by the error maps (next 
section). The results are discussed in Sect.|4] 



2.3. Error Analysis 

Attention has been paid to the error budget in the final images. 
The resultant error is the quadratic sum of individual errors, 
consisting of a) the readout and photon noises added quadrati- 
cally and calculated by the software, denoted Vs,,,/ in Table[5] 
b) an error on the background level which was subtracted; this 
is taken as the l<x dispersion in the histogram of values used 
to determine the sky value to be subtracted and includes un- 
certainties due to residual dark current, flat field, glitches, rem- 
nant images, and long-term drifts, denoted 'l<x Sky Error' in 
Table [3] and c) the effects of the non-stabilization of the sig- 
nal (memory) as a function of time over a temporal frame-by- 
frame block at any pointing; this error has been estimated from 
variations in the signal that are greater than 3cr within a block 
(assuming that such larger variations are caused by only this 
error) and applying a factor which accounts for the fact that the 
error is worse near the beginning of a block (see Galliano 2004 
Roussel et al. 2001al, denoted 'Transient Rms' in Table|5] The 
quadratic sum of these errors is referred to as the random er- 
ror (V in the table) and has been computed 'per pixel', where 
the pixel is the originally sized pixel given in the table, and 
also has been computed 'per square arcsec', after interpolating 
onto a 1 arcsec pixel grid. We formed a total Random Error 
map (showing cr, for each map point, i, from the quadratic sum 
of the <TR ms map, the Transient Rms map, plus a constant lcr 
Sky Error value. These maps are shown in Figs.[Q) through|5]3. 
However, for the beam-switching observations, since we also 
have the sky background errors for each pixel (rather than the 
same sky value for every pixel), all three sources of error were 
calculated in one step (thus the missing intermediate step val- 
ues in Table |3j and the final Random Error map is shown in 
Fig. 05. 

The calibration error (see Table|3j has not been included in 
Figs. \3p through |5j> since this would globally raise or lower 
the flux level of the map and these maps are mainly meant 
to show variations across the fields. However, it is included 
in any calculations of flux density. A remaining error, that of 
variations in the flux density calibration according to orbit (« 
5%, Blommaert et al. 2003 1 has not been included. Since 6.8N, 
7.7N, 9.6N, and 1 1.3N observations were all taken at the same 
orbital point, this error does not enter into a field by field com- 
parison of these 4 bands. However, it will be a factor in a com- 
parison between these bands and the remaining two. 

We note that the true errors cannot be determined rigor- 
ously, for example, the memory effects of the camera are not 
accurately known (Roussel et al. 200 lai, although they can be 
tested for, as described below. Thus the quoted errors, though 
calculated rigorously, should be considered indicative. An ex- 
ample is seen in Fig.[3t in which there remain three residual ar- 
tifacts to the west of the galaxy which correspond to the edges 
of the three fields. We have also, therefore, shown the lowest 
contour level in Figs.^ through^ as twice the rms level on 
the emission map itself, the latter being an indication of the 
pixel-to-pixel variation across a field in the region of the sky. 

Since we are reporting the discovery of halo PAHs in 
NGC 5907 in bands, 7.7N, 1 1.3N, and 6.7W (see Sect.O, we 
performed some additional tests which were designed to check 
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Fig. 1. Maps of NGC 5907 in the 6.8N band showing (a) total emission and (b) the Random error map (see Sect.|2ji. The greyscale 
is in units of /^Jy arcsec~ 2 and ranges from the minimum to maximum map values. Contours are at 0.012 (2a map ), 0.018, 0.03, 
0.06, 0.10, 0.20, and 0.30 mJy arcsec~ 2 where cr map is the rms noise level as measured from the emission map itself. 

Fig. 2. As in Fig.[T] but for 7.7N and with contours at 0.006 (2a map ), 0.01, 0.018, 0.05, 0.10, 0.30 0.50, and 0.65 mJy arcsec" 2 . 
The inset shows the emission after an estimate for the stellar contribution has been subtracted (see Sect. 13. U . 

Fig.3. As in Fig.[T] but for 9.6N with contours at 0.015 (2o- map ), 0.025, 0.035, 0.05, 0.10, 0.18, and 0.25 mJy arcsec 2 . 



whether this high latitude emission might have been produced 
artificially. These tests are as follows, foccusing on the 7.7N 
band. 

We first examined the point spread function (PSF) whose 
FWHM values are 7.2" in these bands (Table 0. The PSFs 
at different locations on the image in these bands have been 
carefully examined and characterized by fitting 2-D gaussians 
(see Galliano 2004 1 and are, to a high accuracy, symmetric. 
Thus, the wings of the PSF should be completely negligible 
7.2" beyond the projection of the disk. To check this, we took 
the 7.7N image (Fig.13, measured its semi-major axis (245"), 
and created an inclined (86.5°) thin disk model (see Irwin et 
al. ^999 for similar modelling). The resulting projected half- 
width perpendicular to the plane at the center of the galaxy is 
15" and less elsewhere. This model was then used as a tem- 
plate to blank all emission outside of the inclined thin disk re- 
gion from the Fig.[2]image. The resulting blanked image then 
contained only emission from the disk region. The blanked im- 
age was then convolved with a 7.2" gaussian PSF to create an 
image consisting of disk plus extended emission due to the PSF 
only. This resulting map was then compared with the original 
image by taking the ratio of the two and determine the pencen- 
tage change per pixel that results. We find that, outside of the 
projection of the thin disk, all ratios are < 1 % and the mean is 
< 0.05%, verifying that the PSF contributes negligibly to the 
halo emission. 

Secondly, we compared results between using the library 
flat field (as described in Sect. 12. 21 and using a sky flat, formed 
from a portion of the image in which the sky is visible. The 
mean relative error between the two results over comparable 
regions is 4.5% and high latitude structure, most of which is 
in the region in which the sky flat could be formed, remains 
present. Moreover, in both images, artifacts, such as those at 
the edge of the field in Fig. [3^, cannot be seen. 

Finally, we checked further into issues related to memory. 
Although the resulting map already has memory corrections 
in place including blanking the first frames of any scan, we 
applied a more aggressive scheme. This involved substantial 
blanking of the first time period when a new scan began, de- 
termining the median of the last 5 frames of the scan, pixel by 
pixel, and then fitting a slope to every pixel independently to 
align the other frames with this median. Some differences in 
structure were noted at a low level, however the high latitude 
emission remained. 



3. The MIR Spectrum and Contributors to the 
Observed Bands 

To interpret the emission in the observing bands, it is first nec- 
essary to understand what the various contributors are to the 
MIR emission. For this purpose, we provide ISOCAM CVF 
(circularly variable filter) spectra of regions in two other galax- 
ies for comparison (see Fig. [7}. The main spectrum is of the 
nucleus of the nearby galaxy, M 83, and the inset is of a quies- 
cent region away from the starburst in M 82. Various features 
are labelled as well as the bands used in our observations. 

There are two main contributors to the MIR spectrum: 
1) strong spectral peaks due to very small grains or large 
molecules (PAHs) which, due to their small sizes (typically 
~ 0.001 /mi) and thus, low heat capacity, are stochastically 
heated, undergoing wide temperature fluctuations covering 
100s of K from absorption of a single photon (Leger & Puget 
IT984I Sellgren et al. [T990l Allamandola et al.Q989l. The pri- 
mary emission bands, relevant to the wavelength range ob- 
served for NGC 5907, are located at 6.2, 7.7, 8.6, 11.3 and 
12.7 /mi and can be seen as strong emission peaks in the spec- 
trum. 2) a continuum composed of very small (of order ~ 0.01 
yum) grains (VSGs) as modeled in our galaxy by Desert et al. 
(1990). Depending on their size and the local radiation field, 
these grains could be stochastically heated or in thermal equi- 
librium. A modeled continuum is shown in the Figure inset, in- 
creasing in strength to longer wavelength. Finally, not shown in 
the figure, is a possible weak stellar contribution that, if present, 
would contribute at the low wavelength end of the spectrum, 
decreasing at higher wavelengths. Since we are able to make 
numerical estimates of this contribution from previous obser- 
vations of NGC 5907 itself, we deal with this separately in 
Sect. El 

Without doubt, the most overwhelming emission in the 
MIR band is due to PAH features. These features are ubiqui- 
tous in normal spirals and account for almost all of the MIR 
energy (e.g. Vogler 2005 Genzel & Cesarsky 2000 & refer- 
ences therein). It is also now established that most normal spi- 
ral galaxies show qualitatively little difference in their spectra 
in the MIR and that the spectral shape is largely independent of 
star formation (Lu et al 2003 1. From higher resolution Galactic 
observations, we know that the VSGs (producing the dust con- 
tinuum) peak in the nebular regions, while the PAHs peak out- 
side the HII regions, in the photodissociation regions (PDRs) 
around molecular cores (e.g. Cesarsky et al 1996, Verstraete et 
al 1996 1. While the PAHs emit profusely in the disks of galax- 
ies, excited predominantly by UV photons, they are also ob- 
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Fig. 4. As in Fig.EI but for 11. 3N with contours at 0.016 (2cr map ), 0.022, 0.038, 0.055, 0.10, 0.20, 0.30, and 0.45 mJy arcsecT 2 . 

Fig.5. As in Fig.[T] but for 12W with contours at 0.028 (2o- map ), 0.05, 0.10, 0.20, and 0.35 mJy arcsec~ 2 . 

Fig. 7. Comparison MIR spectrum of the central position of M 83 from the ISOCAM C VF observations of Vogler et al. ( 2005 1. 
The major features are labelled and the bands used in the observations of NGC 5907 are indicated at the top. Inset: MIR spectrum 
of a quiescent region in M 82, away from the central starburst, showing modeled components, from Laurent et al. ([2000l. 



served in the more diffuse ISM (e.g. Chan et al. 120011 Mattila et 
al. 1996 1 where optical photons are thought to be the source of 
excitation. Optical photons are also an explanation for the pres- 
ence of PAHs in elliptical galaxies (Athey et al 2002 Xilouris 
et al 2004 1, where evolved stars are the primary stellar popu- 
lation. Moreover, in a study of 5 galaxies over a range of star 
forming activity, Haas et al. (2002i find a spatial correlation 
between PAHs and the cold dust distribution which is more 
widespread than around SFRs alone. 

We thus expect the relative importance of the continuum 
to diminish in galaxies of lower star formation rate or in qui- 
escent regions of galaxies. This is indeed observed in M 83 
(the increasing continuum longward of A14 fim at the nucleus 
shown in Fig.[7]becomes much less prominent in the interarm 
regions) as well as other known quiescent star forming regions 
(Fig. □ inset). Thorough modeling of ISOCAM CVF spectra 
using Lorentizian profiles further supports the dominance of 
the PAH bands, greatly reducing the need for a strong contin- 
uum. (Boulanger et al. 1998 1. Indeed the three PAH features 
at A6.2, 7.7, and 8.6 fim in M 82 can be fit without any sig- 
nificant underlying continuum at all; a continuum contribution 
is fit with a simple linear curve, increasing with wavelength 
(Laurent et al. 2000 1. Vogler et al. <2005> also find that a con- 
tinuum contributes only ~ 5% for these three PAH bands. As 
this is less than a typical error bar on the observations, we con- 
sider the continuum contribution below these PAH features to 
be negligible. 

Similar arguments apply to the ionic emission lines. In nu- 
clear regions, a [Nell] emission line at 12.8 fim has been ob- 
served and is blended with the 12.7 fim PAH feature. However, 
in a sample of 69 nearby normal spiral galaxies with ISOCAM 
spectra, Roussel et al. ( 2001c I found that the contribution from 
the [Nell] line is negligible. The A7.Q fim [Aril] line, visible at 
the nucleus of M 83 (Fig.[7J is no longer seen in the interarm 
regions of this galaxy (Vogler et al. 12005b . 

The final possible MIR feature to be remarked upon is a 
broad silicate absorption feature at 9 fim which has been de- 
tected in narrow central regions of very dusty starbursts and 
AGNs (e.g. Laurent et al. 120001 Tran et al. 12(31)11 Sturm et al. 
2000 1. This feature is not likely to be important toward the 
central regions of normal starbursts or in quiescent disks (e.g. 
Roussel et al l2001all200Tbl Vogler et al l20031 Dale et al. l2UoT1 
also arguments in Spoon et al. 2004 1. even in the case of edge- 
on spiral galaxies (Mattila et al. 1 19991 for NGC891). 

Therefore, the MIR spectrum of a quiescent galaxy like 
NGC 5907 can be expected to be dominated by the PAH emis- 
sion bands, with a weak underlying continuum which becomes 



stronger at longer wavelengths. We thus interpret the emission 

in the various observed bands as follows: 

6.8N: Lorentzian wings of the A6.2 fim and A7.7 fim PAH 

bands. 

7.7N: A7.7 fim PAH band. 

9.6N: Lorentzian wings of the A7 .7 fim, /i8.6 fim, and Al 1.3 fim 
PAH bands plus possible minor contribution from a continuum. 
11.3N: Al 1.3 fim PAH band plus possible continuum. 
6.7W: Combined A6.2, A7.7, and part of the AS.6 fim PAH band 
emission. 

12W: Broad band emission, equivalent to the IRAS A 12 fim 
band, dominated by PAH features with likely continuum. 

3.1. The Stellar Contribution 

The stellar contribution to the MIR band is illustrated beauti- 
fully in Lu et al. d2003t for both a reddened and de-reddened 
case (their Fig. 6). These authors have obtained complete ISO 
spectra in the 2.4 to 5.9 fim and 5.8 to 11.6 fim bands for 45 
galaxies finding, as have other authors, that the characteristic 
shape of the spectrum in these bands is extremely consistent, 
galaxy to galaxy. The stellar contribution can be approximated 
by a modified black body with a A~ 2 emissivity law. This be- 
haviour for the stellar contribution also accurately describes the 
elliptical galaxies in their sample. Thus, we expect the stellar 
contribution in both the disk and halo regions of NGC 5907 to 
have a similar A dependence. 

We have obtained the K, band (A 2.159 fim) image for 
NGC 5907, from the 2MASS Large Galaxy Atlas (LGA, Jarrett 
et al. 2003 1, which is the band most likely to show only stars 
with negligible dust obscuration. This image was interpolated 
onto the same grid as our MIR images and then smoothed to 
7.2" resolution. The 'data numbers' (dn), in the FITS file were 
converted to Jy using the transformation, f v = 6.248 x 10~ 6 dn 
Jy 2 . As a check on our processing steps, we measured the to- 
tal flux density of the source in the regridded, smoothed image, 
finding f v = 1.322 Jy which corresponds to niK s = 6.76. This 
agrees with the magnitude given in the LGA for NGC 5907 of 
m.K s = 6.76 + 0.02. The resulting image has an rms noise of 
o" = 1.6 fiJy arcsec~ 2 . From the emissivity law (following Lu 
et al.), we expect the stellar contribution, extrapolated into the 
MIR bands in which PAHs dominate to be 7.9% in the 7.7N 
band, 3.7% at 11. 3N, and 10% at 6.7W of whatever value is 

2 This transformation uses the calibration, m Ks = ZP - 2.5 log(dn) 
and f v = / v (0) lfr 04 "'*' where ZP = 20.0704994 is the zero point 
provided in the FITS header, and /,,(0) = 666.7 Jy, the latter from 
http://ssc. spitzer. caltech. edu/documents/cookbook/html/cookbook. html 



8 



Judith A. Irwin and Suzanne C. Madden: Discovery of PAHs in the Halo of NGC 5907 



in the K s band at the relevant location. With this extrapolation, 
we created 3 maps from the K s band image for the 3 bands 
dominated by PAH emission and subtracted each from the total 
emission map. These results are shown in the insets to Figs. [3] 
0] and [6] with the same contour levels as the main maps. We 
then measured the global flux of the stellar-subtracted maps 
and compared this to the global flux of the main maps finding 
a difference of 5%, 2%, and 2% for the 7.7 /jm, 1 1 fim, and 6.7 
pm maps, respectively (ignoring emission from the two outer 
frames in the latter case). Since the flux corrections are well 
within other error bars (Sect.|5J and the appearance of the maps 
is also close to the original non-subtracted maps, we consider 
the stellar contribution to be negligible and continue our dis- 
cussion and analysis with respect to the unsubtracted maps. 

4. Results 

All maps along with the best estimated error maps are shown 
in Figs.[0to[6]as described in Sect. 12.21 The error maps delin- 
eate the final field of view. Emission has been detected in every 
observing band as shown in these figures. 

4.1. Disk Emission 

Emission is clearly detected along the disk of NGC 5907. In no 
individual band do we have full coverage of the entire galaxy. 
However, for the narrow bands, 6.8N, 7.7N, 9.6N, and 11. 3N, 
we have full coverage of the south-east major axis and for 
6.7W, small segments at both ends of the galaxy have been ob- 
served. In all cases, the emission extends along the optical disk 
with detections to the 'end' of the optical disk in some cases. 

The extent of the PAH emission, specifically, is illustrated 
by Fig. [S^ in which we show the smoothed 7.7N emission (A 
7.7 fim PAHs) over the DSS optical image. It is clear from this 
overlay that the PAH emission extends all along the optical disk 
of the galaxy. It is further illustrated in Fig.[9]in which we show 
the smoothed, high sensitivity 6.7W emission (A6.2, 7.7, and 
8.6 fim PAHs) superimposed on the HI total intensity map from 
Shang et al. ( 1998 1. PAH emission is seen very far out along the 
galaxy disk where the HI emission is also strong. The total ex- 
tent of the PAH emission along the major axis is 11.0 arcmin 
(~ 35 kpc), as measured to the outermost detectable contigu- 
ous features on the unsmoothed map (cf. 41 kpc for the optical 
disk, Table [0. However, there is also some evidence for real 
features farther out. For example, several emission peaks can 
be seen at the edges of the 6.7W field of view on the south-east 
end of the major axis at RA = 15 h 16 m 12 s , DEC = 56° 13' 
09". A channel-by-channel comparison (in velocity) of the HI 
distribution with the PAH emission, however, fails to show any 
clear correlation. 

The lack of global correlation with HI is further illustrated 
in Fig. 110b in which we show major axis slices of the HI to- 
tal intensity distribution, the CO, and the smoothed 7.7N emis- 
sion. It is clear from this plot that the PAH emission follows the 
molecular gas distribution and not the atomic gas distribution. 
This was also found by Dumke et al. ( 1997 1 for the Al. 2 mm 
dust distribution as well as by Alton et al. (I2004> for the A850 
fim emission. Fig.lToh shows two sets of peaks in the A7.7 fjm 



Fig. 8. (a) The 7.7N band image (A7.7 fim PAH band), 
smoothed to 12" resolution superimposed on the optical DSS 
image. Contours are at 4.6 (2o- map ), 7.5, 13, 25, 50, 75, 200, and 
380 x 10~ 3 mJy arcsec~ 2 . (b) The 7.7N band image as in (a) su- 
perimposed on a greyscale /1850 SCUBA //m map from Alton 
et al. d2004> with a resolution of 16". Note that the SCUBA 
field of view ends at a Declination of 56° 17'. 

Fig. 9. The 6.7W band image, smoothed to 12" resolution su- 
perimposed on the HI total intensity image (data from Shang 
et ai. [T998l Contours are at 0.0004 (2cr mflp ), 0.0008, 0.0015, 
0.003, 0.01, 0.04, and 0.10 mJy arcsec 2 . 

distribution on either side of the nucleus, one at ~ + 70" and 
the second at ~ + 100". The CO comparison also show two 
peaks, interpreted by Dumke et al. ( 1997 1 to represent rings or 
spiral arms, but at somewhat offset radii of + 60" and ± 120". 

Fig. 10. Total intensity cuts along the major axis of NGC 5907. 
Positional uncertainties are of order +7". (a) The HI total inten- 
sity (faint extended curve) from data of Shang et al. ( 1998), the 
CO distribution (faint solid plus dashed curves showing data 
points and error bars) from Dumke et al. ( 1997 1, and the Al .1 
/urn from this work (darkest curve), the latter smoothed to the 
same resolution as the HI data (18.6" x 17.8" at PA = -38°). 
The CO resolution is 21" and the center point has been placed 
approximately at the CO center of RA(J2000) = 15 h 15 m 53.8 s , 
DEC(J2000) = 56° 19' 37". (b) The ^7.7 /jm PAH emission 
from the 7.7N band, smoothed to 16" resolution shown with 
the /1850 //m SCUBA emission at the same resolution. Curves 
and intensities represent an average over a 16" swath parallel to 
the major axis. The arrow marks the edge of the SCUBA field 
of view. 

In Fig.llOb. we show a comparison between major axis cuts 
for our ISO A1J fim data and the /1850 //m SCUBA data, both 
averaged over a minor axis extent of 16" and at the same spatial 
resolution. The shape of the PAH emission curve and the /1850 
//m emission curve, the latter which traces cool dust, are re- 
markably similar showing the same number and approximate, 
but not exact, positions of peaks on either side of the nucleus. 
Positional offsets between between the two different slices for 
5 peaks range from zero to 19" with a typical value of ~ 13". 
Given that the resolution is 16" and relative positional uncer- 
tainties between maps is of order 7" (see caption), the positions 
of the peaks are likely consistent, though we cannot rule out 
the possibility of some displacement which might be detected 
with higher resolution observations. The line ratios from these 
curves (calculated though not plotted) are also consistent with 
the values of Haas et al. ( 2002 1 who found A7.1 (im/A850 fim = 
2 (1.8 to 2.2 for 5 galaxies). Their denominator is from SCUBA 
observations but the numerator represents the peak A7.7 fim 
flux rather than the average in the 7.7N band as we have plot- 
ted. All of the emission in our 7.7N band is from the Al .1 fim 
PAH line but the peak, if we could measure it directly, would be 
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higher than the average. From Fig.0 the increase should be of 
order a factor of a 1.5 which would put our ratio in even closer 
agreement with Haas et al. Therefore, overall, our results are 
consistent with a possible relationship between PAHs and the 
850 fim cool dust distribution. 

Finally, there is excess A7.7 fim emission between radii of 
-150" and -250" (Fig.llOh) in comparison to the CO distribu- 
tion. (Note that we cannot discern this in Fig.llOb because of 
the truncated field of view of the SCUBA map.) We believe 
that this excess is significant. Indeed, Dumke et al. ( 1997 1 also 
found an excess in the A1.2 mm emission in comparison to CO 
over the same radii but on both the north and south sides of the 
nucleus. This shows that there is PAH emission as well as cool 
dust where there is little molecular gas. We will return to this 
point in Sect. [5] 

In Fig.^2 we show the broad band 12W emission (equiva- 
lent to the IRAS 12 fim band) superimposed on the DSS opti- 
cal image. This MIR band should contain the highest contribu- 
tion of continuum emission from VSGs compared to the other 
bands (Fig.0 inset). The 12W band has the highest spatial res- 
olution of all observations (1") and reveals that the strongest 
ridge of MIR emission lies along the optical dust lane (note 
that the field of view does not cover the entire disk), indicating 
a good correlation between VSGs and the optically obscuring 
dust. 



Fig. 11. The IRAS band-equivalent image, 12W, (contours) 
over the Digitized Sky Survey (DSS) optical image. Contours 
are at 0.028, 0.050, 0.10, 0.18, 0.28, and 0.36 mJy arcsec~ 2 . 
The 12W field of view does not extend over the whole optical 
image (see Fig. |5j>). The known faint stellar ring (Sect. Il.lt . 
which cannot be discerned in this image, intersects the disk at 
declinations of 56° 23.4' and 56° 16.8' in projection (Zheng et 
al. rP?99l . 



4.2. High Latitude Emission 

A particularly striking result is the evidence for high latitude 
emission, especially in the 7.7N band (Fig.|2), the 1 1.3N band 
(Fig. @}, and the high-sensitivity beam-switched 6.7W band 
(Fig-EJ which all trace PAH emission. The 1 1.3N band in prin- 
ciple, can contain PAHs plus continuum, but the continnuum 
should be negligible so far from star forming regions. Thus it 
is clear that PAHs exist far from the plane of NGC 5907, in the 
case of 6.7W, as far as 6.5 kpc (2') from the plane, though ~ 3 
kpc is more typical. Even the 6.8N band (Fig.[0, which should 
consist only of emission from PAH band wings, shows some 
evidence for features away from the plane. Given the galaxy's 
edge-on orientation (86.5°, Tabled, the correction for galaxy 
inclination is negligible. For example, for a PAH disk radius of 
250" f Fig. llOb . the semi-minor axis is only 15.3" in a thin disk 
model. Moreover, the CO semi-minor axis is only ~ 4 arcsec if 
we use the HWHM of 200 pc modeled by Dumke et al. Jl997t . 
These features resemble the extra-planar emission seen in other 
galaxies, (see, e.g. the radio continuum features in NGC 5775, 



Lee et al. 2001 1 and some are similarly arc-like in appearance 
(e.g. the east side of the A7.7 yum emission in Fig.|2j. Recall also 
(Sect. rTTl that the stellar scale height has been measured to be 
0.49 kpc (9"). Note that the narrow bands, 7.7N and 11. 3N, 
show more structure than the wide band, 6.7W which may re- 
sult from the ability of these bands to isolate a single PAH line. 
Although structural details may vary, both narrow bands show 
high latitude features at approximately the same locations (see 
also Sect.l2"3l. 

Given that the beam-switched 6.7W band observations are 
our most sensitive and also the fact that this band is dominated 
by several PAH emission features, we have taken the origi- 
nal, unsmoothed 6.7W data shown in the central regions of 
Fig- ED and averaged the emission in strips parallel to the ma- 
jor axis. The resulting minor axis profile is shown in Fig. 1121 
The wings on this averaged profile approach the noise between 
100" and 150" (5.3 kpc to 8.0 kpc) and the line half-width at 
the 3cr level is ~ 95" (5 kpc). We fit various functional forms 
to the profile up to 3 components in total, including Gaussians, 
Lorentzians, Voigt, combinations of these forms and a com- 
bination of Gaussian + exponential. The best fit results were 
achieved by 3 Gaussians or 1 Gaussian plus 1 Voigt profile, 
which resulted in equivalent residual rms values. However, 
since the Voigt profile does not have easily identifiable param- 
eters, we show the fits (Fig. ll2t for the Gaussian combinations 
only and provide tabular information (Table 0J| for these com- 
binations as well as the Gaussian + exponential for compara- 
tive purposes. The broadest Gaussian scale length is 65.3" (3.5 
kpc). Increasing the number of components would reduce the 
residual rms further but such a situation is likely not physical. 
The profile, in general, consists of a central Gaussian compo- 
nent plus broad wings. Also, since the high latitude emission is 
highly structured (Fig. |6j, the scale length of the broad com- 
ponent will vary with position. Nevertheless, it is clear that 
an extended halo or thick disk of PAH emission is present in 
NGC 5907 with a characteristic scale height of 3.5 to 5 kpc. 

Fig. 12. Average profile of the PAH emission as a function of 
z from the high sensitivity 6.7W data (Fig.|6|l. Negative values 
of z represent the western side of the major axis, (a) Plot of 
the data over the full range of intensity, (b) Plot of the data 
and two of the models over only the lowest 1.7% of the total 
intensity. The dark solid curve shows the data with the noise, 
at an rms level of 2 x 10~ 4 mJy arcsec -2 , now obvious. The 
grey curve shows the best fit 2-component gaussian model and 
the narrower dark solid curves shows the best fit 3-component 
gaussian model. See Table@]for the parameters of the gaussian 
fits. 

Given the apparent correlation between the ISO A7.7 yum 
emission and the SCUBA /1850 yum emission in the disk 
(Sect. 14. H . it is of interest to see if this correlation extends 
to the halo region. Fig. [SJ> shows a comparison between these 
bands. Note that the /1850 fim field of view is truncated at a 
Declination of 56° 17' and the 7.7N map is truncated on the 
eastern side. Note also that the greyscale showing apparent 
emission away from the plane in the SCUBA data is all below 
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Table 4. 6.7W z Profile 



Data or Type of Fit 


Component No. 


Peak Intensity 


FWHM a 


Residual rms 






(mjy arcsec~ 2 ) 


(arcsec) 


(mjy arcsec -2 ) 


Data 








0.20 x 10" 3 


1 -component Gaussian 


1 


0.315 


18.3 


3.96 x 10- 3 


2-component Gaussian 


1 


0.287 


16.2 






2 


0.0376 


43.2 


1.07 x 10~ 3 


Gaussian + Exponential 


1 


0.264 


16.6 






2 


0.065 


16.5 b 


1.11 x 10~ 3 


3-component Gaussian 


1 


0.257 


15.5 






2 


0.0564 


27.2 






3 


0.00131 


65.3 


0.86 x 10- 3 



a Unless otherwise indicated. 
b Exponential scale length. 



the 3cr level (Alton et al. 2004 their Fig. 1) with the exception 
of the 'detached' feature at RA = 15 15 52, DEC = 56 17 45. 
Indeed, Alton et al. find an exponential scaleheight of only 0.11 
kpc (2") for the /1850 //m map. If the AE50 /im emission existed 
in the halo at the same fraction of the in-disk emission as the 
PAHs, then the A850 fim emission should have been detected in 
the halo above the noise of the Alton et al. map. This suggests 
that either the large grains are under-represented in the halo in 
comparison to the PAHs, or are under-emitting. 

The broad band 12W emission shows very little vertically 
extended emission (Fig. [5] Fig. llH in comparison to the PAH 
bands (cf. Figs.|5]and|3. This is likely due to the fact that the 
dynamic range in the 12W band is only 10 to 17% as great as in 
the 7.7N and 1 1 .3N bands (Tabled. For example, the extended 
vertical emission seen in Figs. 13 and |4] present in Fig. [5] at 
the same relative brightness with respect to the maximum disk 
emission, would mostly fall at the level of the noise or lower. 
One obvious feature is visible at the edge of the field of view 
on the north, extending to the east. (The optical feature within 
this extension appears to be a foreground star.) The reality of 
the large scale feature is not clear but its 'footprints' close to 
the disk have counterparts in the PAH bands (e.g. Fig.0} sug- 
gesting that the PAH contribution within this band may be re- 
sponsible for it. 

4.3. Flux Densities and Band Ratios 

In no case do we have complete coverage of the galaxy so 
global flux densities cannot be obtained. However, we have de- 
termined the flux density for each of the fields displayed. These 
are listed in Table|5] Since the 6.8N through 1 1.3N bands have 
the same field of view, their fluxes can be compared directly. 
There is one common field for all maps, however, which is out- 
lined in Fig.|6]3. Thus, we also calculate fluxes for this common 
region (Table |5} and compare these results for all bands. The 
resulting flux ratios are provided in Table|6]with some compar- 
ative values, where known, also given. 

The 6.7W:12W ratio is interesting since it is significantly 
higher than the values given by Dale et al. J2001i . This ratio 
represents a comparison between a broad ISOCAM band that 
contains mainly 3 PAH bands (Sect.|5J with the broader 12 jjm 



IRAS -equivalent band which may contain both PAHs and con- 
tinuum. Dale et al. (2001 1 determined average values of this ra- 
tio, using the IRAS 12/im flux itself, for normal galaxies in dif- 
ferent /eo/Zioo bins. The IRAS flux ratio, /60//100 is a measure 
of the heating of classical grains in temperature equilibrium. As 
such, it is widely used as a measure of star formation activity, 
with an increasing ratio implying greater dust heating (Helou 
IT9861 Dale et al. l2000ll200lt . As the SFR increases (/eo//ioo 
increases), the 6.7W/(broad-band 12jt/m) ratio decreases, re- 
flecting the fact that a) the contribution of a hot dust continuum 
increases with SFR and b) PAHs are expected to be destroyed 
in regions of high heating intensity such as active star forming 
regions (Dale et al. 2001 and references therein). The /eo/Zioo 
ratio for NGC 5907 (0.24, TableQ places it at a SFR that falls 
below all of the galaxies of the Dale et al. sample. Moreover, 
the /60//100 ratio for NGC 5907 could be even lower than that 
quoted in Tabled i.e. 0.21 from earlier ISO data, Bendo et al. 
2002 Thus, the high 6.7W: 12W ratio for NGC 5907 (Tableg} 
is likely consistent with the data from other galaxies in that it 
would represent an extension to lower SFR. 

The ratios for the bands 6.8N through 11. 3N (TablegJ all 
have the same field of view and so can be compared directly. 
(Note also, that all results for the 'field' and for the 'com- 
mon' regions agree within errors.) For these bands, we make 
a comparison to the galaxy, M 83, for which similar obser- 
vations have been made. NGC 5907 is not perfectly identical 
to M 83 in its global properties (Table The two galaxies 
are of similar morphological types and have similar global star 
formation rates, infrared luminosities, and modelled tempera- 
tures. However, M 83 is physically smaller than NGC 5907. 
Thus, NGC 5907 contains an order of magnitude more dust but 
M 83 is the more active galaxy when SFR per unit area is de- 
termined. Nevertheless, the shapes of the different spectra in 
M 83 (see Vogler et al. 2005 and Fig. [7) are remarkably con- 
stant as a function of location in that galaxy, i.e. as a function of 
differing local SFR, the main differences being increased con- 
tributions of continuum, especially longwards of 14/im (which 
we sample only marginally in NGC 5907 in the 12W band). 
For M 83, the main differences between spectra can be seen as 
an increased continuum contribution from hot grains which be- 
come more prominent towards the nucleus. Thus, we provide 
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Table 5. Flux Densities 



Band 


Flux Density (Field)" 


Flux Density (Common)" 




(mly) 


(mJy) 


6.8N 


1170 ± 97 


519 ±42 


7.7N 


2170 ± 141 


918 ±60 


9.6N 


919 c ± 66 


360 ± 26 


11.3N 


2071 ± 116 


828 ± 46 


12W 


940 ± 38 


557 ± 57 


6.7W 


666 ± 22 


591 ± 20 



" Flux density determined over the entire field shown in Fig. through Fig.|6|i. The error (here and throughout) is calculated via y£of 
over the map pixels, i, where cr, is the quadratic sum at each pixel of a rms , the Transient rms, the ltr sky error and the calibration error (see 
Table^}. 

b Flux density determined over the region that is common to all maps and outlined in Fig.|6j). 
c Excluding the 3 small regions with residual errors (see Sect. l2.3> . 



Table 6. Flux Density Ratios 



Bands" 



Ratio or Relation 



[0.52!^]: 1" 
[1.06±0.16] c :l 
[0.57 ± 0.08]:1:[0.39 ± 0.05]:[0.90 ±0.11] 
[0.54 ± 0.08]:1:[0.42 ± 0.06]: [0.95 ± 0.12] 
0.49 : 1 : 0.59 : 1.12" 
0.55 : 1 : 0.48 : 0.95 e 



Galaxy or Region 



Ref. 



6.7W:12(IRAS) 
6.7W:12W 

6.8N:7.7N:9.6N:11.3N 
6.8N:7.7N:9.6N:11.3N 
6.8N:7.7N:9.6N:11.3N 
6.8N:7.7N:9.6N:11.3N 



8 galaxies 
Common d 
Common d 
Field d 
M 83 (R < l') f 



Dale etal. 2001 
This Work 
This Work 
This Work 
Vogler et al. <2005l 



M 83 (1' < R < 3') f Voeler et al. <2005l 



" The center wavelengths of the bands is given below each sequence that is the same. 

b Average of 8 galaxies in a flux ratio, /60//100, bin from 0.28 to 0.35 (Dale et al. 2001 1. 

c This error includes an additional 5% error to account for variations in orbit (Sect. l2"3l . 

d Common or Field region as quoted in Table|5| 

e Errors (estimates only) are believed to be of order 20%. 

f R is the projected radial distance from the center of the galaxy. 



the 6.8N through 11. 3N ratios for M 83 for the inner, hotter, 
1' radius region as well as the outer, cooler, region between 
radii of 1' and 3' in Table [6] Given the estimated errors on the 
various ratios, there is agreement between all values for M 83 
and those of NGC 5907, although the most variation is seen for 
9.6N. The band, 9.6N, measures a mixture of PAH wings plus 
possible hot continuum. The 9.6N:7.7N ratio for NGC 5907 is 
somewhat low in comparison to the nucleus of M 83 but agrees 
with the result for the interarm region of M 83. This is again 
quite consistent with the low SFR for NGC 5907. 

The Dale et al. comparison values for 6.7W:12W are glob- 
ally determined values. Our data allow us to further examine 
this ratio as a function of projected distance from the major 
axis, as shown in Fig.^] In Fig. 113b. we show the increasing 
6.7W:12W ratio with projected distance from the major axis. 
Given that this ratio can only be determined over a restricted 
area of the 'common' region due to signal-to-noise consider- 
ations, this may reflect either an increase in this ratio with z 
height or an increase with distance from the center of the galaxy 
since the region shown does not extend beyond the projected 
optical disk. In either case, this is likely due to an increase in 
the distance from star forming regions and therefore a decrease 
in the contribution of the hot continuum that would be expected 
in the 12W band. In Fig. 113b. we show the ratio of the emis- 



Fig. 13. (a) Profile of the 6.7W/12W ratio from the common 
region only (see Fig. |6p) as a function of projected distance 
from the major axis. The values are averages over 71" strips 
parallel to the major axis, (b) The z profile of the 11.3N/7.7N 
ratio for the fields shown in Figs. [2] and |4j but averaged over 
320" strips. Note that the z extent is much larger than in (a). 



sion in two PAH bands: 1 1.3N/7.7N. In this case, the ratio can 
be determined to much higher latitude. The maximum change 
in ratio is of order 2.4 with the in-disk value lower. Given that 
more continuum, if present, is expected in the 1 1.3N band than 
the 7.7N band in the disk, then the in-disk ratio could be lower 
than shown. As indicated above, the global ratio of these two 
bands is quite consistent with what is seen in other galaxies so 
we may be observing a real variation in the 11.3/7.7 fj.m PAH 
ratio between the changing environments of the disk and halo. 
The nature of this variation is unclear and the exact magnitude 
will require the acquisition of spectra and modelling. However, 
we note that the variation seen here is not inconsistent with the 
variations between PAH feature strengths observed for regions 
in our own Milky Way, i.e. of order a factor of 5 (see Peeters et 
al. l2"UU2l 
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4.4. Radio Continuum Correlation? 

Vogler et al. ( 2005 ) find a good spatial correlation between the 
A6 cm radio continuum emission and the 6.7W (largely PAH) 
MIR emission in M 83, and Brar et al. ( 2003 1 also find a good 
correlation between the A850 fim emission and the 617 MHz 
radio continuum emission in the edge-on galaxy, NGC 5775. 
These results, which link the MIR or sub-mm with the radio 
continuum, are reminiscent of the well-known FIR-radio cor- 
relation whose origin is not fully understood but is thought 
to relate to star formation in galaxies (see Groves et al. 2003 
and references therein). Thus, a comparison of our ISO re- 
sults with radio data is of interest. The best available radio data 
sets are those of the NRAO VLA Sky Survey (NVSS) at 1420 
MHz (Condon et al. 119981 and the Westerbork Northern Sky 
Survey (WENSS) at 327 MHz (Rengelink et al. [T997t . The ra- 
dio emission in the in the VLA Faint Images of the Radio Sky 
at Twenty-Centimeters survey (FIRST) is too faint for com- 
parison (but see note below). In Fig. [2] we show the 7.7N 
(/17.7 fim PAH) emission, which is unlikely to contain any sig- 
nificant continuum component, in comparison to these two ra- 
dio images. The strong radio source at the end of the south- 
eastern major axis (actually a double) is believed to be from a 
background source (Dumke et al. 2000). There are significant 
differences between the fields of view and resolutions of the 
ISO and radio data sets which make comparison difficult, and 
smoothing the 7.7N emission results in confusion and trunca- 
tions near the field edges. In addition, details in the structure 
of the 7.7N map at low emission levels may be approximate. 
Thus any apparent correlations are tentative and are presented 
so that higher resolution data in the future might be compared 
to these results. In particular, a radio continuum extension on 
the west side of the major axis at RA * 15h 15m 45s, DEC * 
56° 20' shows extended PAH emission (though the ISO field of 
view also truncates here). South of this feature is another at RA 
=s 15h 15m 50s DEC ~ 56° 18'. On the east side of the major 
axis, there is a feature at RA * 15h 16m 00s DEC ~ 56° 20'. 
These features in the radio are more easily seen in the 327 MHz 
image and at least the northernmost extension on the west side 
can also be seen in a 20 cm image of Dumke et al. (.2000). We 
are not able to determine correlations along the disk due to the 
limited ISO field of view and resolution and also because of 
confusion with the strong background source. 

Finally, we note that there is a previously unreported point 
source in the FIRST data with coordinates, RA(J2000) = 15 
15 49.10, DEC(J2000) = 56 20 27.1, which is displaced 38" 
to the west and 43" to the north of the nucleus (2.0 and 2.3 
kpc, respectively, in projection) with a flux density of 3.1 mJy. 
This flux density corresponds to a spectral power of 4.5 x 10 19 
Watts/Hz which is intermediate between that of Cas A and the 
brightest supernova remnant in M 82 (see Irwin et al. 2000 ) and 
is therefore likely a supernova remnant or remnants. 

5. Discussion 

The mid-IR emission in NGC 5907 is consistent with 
NGC 5907 as a low SFR galaxy. As described in Sect.gj] this 
galaxy is 'PAH-rich' and 'hot-dust-poor'. Our flux ratios are 



consistent with those measured in other galaxies at the low end 
of the range of SFR and are also consistent with those of Vogler 
et al. ( 2005 ) for M 83 for the regions between spiral arms. For 
M 83, the continuum contributes no more than « 10% of the to- 
tal MIR 4 to 18 yum emission with a 90% due to PAHs in these 
regions. Although we do not have CVF data for NGC 5907, all 
other indicators suggest that the situation is likely comparable. 

It is quite clear that the PAH distribution follows that of 
the CO (Pig. flOl, Both the ,1850 fim and A1.2 mm cool dust 
also follow the CO distribution, and the PAH and /1850 fim are 
spatially well correlated within the disk of NGC 5907. This 
suggests that both the PAHs and larger grains are associated 
with the molecular gas component. PAHs are well known to 
trace classical photodissociation regions (PDRs) of molecular 
clouds which are excited by far-UV photons from a hot stel- 
lar population in our own Galaxy (see Sect. [3. However there 
is also evidence that the PAHs are extended radially in com- 
parison to CO both from our results as well as the cool 1.2 
mm dust distribution of Dumke et al. ( 1997 ). The cool dust and 
PAH emission both extend to radii of ~ 250" in comparison to 
~ 175" for CO. While there may be some weak CO present, 
the PAH and cool dust distributions show an excess with re- 
spect to the CO at large radii in comparison to smaller radii. 
These results suggest that, while PAHs are strongly associated 
with the molecular gas, a more widespread distribution, likely 
associated with cool HI, is also present but at a much weaker 
level of emission. 

Haas et al. ( 2002 ) also find a correlation between PAHs and 
cold dust, the latter as measured by the SCUBA /1850 fim con- 
tinuum emission, in galaxies. These authors suggest that PAH 
carriers are widespread and correlate with large dust grains and 
neutral molecular clouds, but require a minimal radiation field 
in order to be excited, thus explaining global correlations be- 
tween SFR and PAH emission as well as the fact that PAHs are 
distributed more widely than the SF regions alone. This is also 
consistent with Peeters et al. (|2004) who suggest that PAHs 
may be better tracers of B stars, which dominate the Galactic 
stellar energy budget, rather than than as a tracer of massive 
star formation. As indicated in Sect. [3] the PAHs in our Galaxy 
peak outside of HII regions but are also profuse emitters in 
the diffuse ISM of galaxy disks, where the excitation may be 
from optical photons from an older population. These studies 
therefore support the association of PAHs with more extended 
neutral atomic gas. Tacconi-Garman et al. (2005) find that the 
global distribution of PAHs in NGC 253 and NGC 1808 fol- 
lows the starburst activity, but on small scales, there is little cor- 
respondence (correlation or anti-correlation) with star forming 
regions. It may well be that PAHs are widespread and excited 
by both UV photons from a hot stellar population as well as op- 
tical photons from the more widespread stellar population from 
which the interstellar radiation field (ISRF) originates. The fact 
that they correlate globally with the CO distribution would then 
be a statement that there are more UV photons available at the 
appropriate wavelength and also that there is a greater concen- 
tration of dust and PAHs in regions of higher density molecular 
clouds. In regions (or galaxies) with very high SFRs the PAHs 
will be destroyed in comparison to the VSGs; this is supported 
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Fig. 14. (a) The 7.7N band image, with contours and field of view as shown in Fig. [2] superimposed on the NVSS A 20 cm radio 
continuum image with low level emission enhanced. The greyscale ranges from 0.01 (« lcr) to 4 mJy beam -1 , with a beam 45" 
in size. The strong source at the south-east end of the major axis is believed to be a background source (Dumke et al. 2000). (b) 
WENSS 327 MHz image (contours) superimposed on a greyscale 7.7N band image. Contours are at 4.5 (1.5<x), 6, 9, 12, 15, 20, 
30, 40, 70, 110, and 160 mJy beam -1 with a beam size of 54 in Right Ascension and 64.9 in Declination. The greyscale ranges 
from 0.003 to 0.02 mJy arcsec -2 . 



by observations of M 83 (Vogler et al. 2005) and of starbursts 
and AGNS (e.g. Sturm et al. l2000t . 

An important result is the presence of PAH emission to very 
high latitude, with scale heights of order 3.5 to 5 kpc and up 
to 6.5 kpc in some locations. By comparison, the CO FWHM 
in the z direction is only ~ 8" = 426 pc after deconvolution 
with the disk (Dumke et al. 119971 and the /J850 //m exponen- 
tial scale height is only 110 pc to the limits of those data. An 
equally significant result is that such high latitude gas is seen in 
a galaxy with a low SFR. To our knowledge, this is the first evi- 
dence for high latitude PAH emission in any external galaxy. A 
A3. 3 //m PAH spur of z extent < 120 pc has now been detected 
by Tacconi-Garman et al. ( 2005 1 in the outflowing wind of the 
starburst galaxy, NGC 253, lending support to the concept of 
'mass-loaded' nuclear starburst winds and the possibility that 
some PAHs can survive such blasts, the intense UV radiation 
and shocks that are expected to accompany them. 

The situation is quite different in NGC 5907, however, 
which does not show strong starbursting or outflowing winds. 
There is, however, a feature with anomalous velocities in high 
resolution CO data that cannot be explained by a bar or rota- 
tion model (Garcia-Burillo et al. 1997 their Feature F). A CO 
loop-like feature of size ~ 20" associated with the anomalous 
structure in the CO maps is suggestive of a small-scale wind 
outflow. The radio point-source described in Sect. 14.41 (not at 
the same location as the CO feature) indicates that some SF ac- 
tivity is indeed occurring now in NGC 5907, although at a low 
level in comparison to other SF galaxies. The main indepen- 
dent evidence for high latitude emission in NGC 5907 is from 
the radio continuum (Dumke et al. 2000 and Fig. I14t which 
shows a A20 cm thick disk scale length of 1 .5 kpc and a large 
scale field of ~ 1 fiG. This indicates that cosmic rays and mag- 
netic fields do exist in the halo, although a detailed search for a 
correlation has not yet been possible (Sect. l4~4l >. 

A mechanism for depositing PAHs in the halo of a low SFR 
galaxy is not yet clear. Since SF and magnetic fields are both 
clearly present in NGC 5907, the PAHs may be vented into 
the halo initially from SNe and stellar winds from SF regions 
and then assisted in reaching high latitudes via coupling to the 
magnetic field. Excitation could be from UV photons that leak 
from the disk into the halo. An interesting possibility is that the 
3-D VSGs are themselves converted into planar PAH macro- 
molecules (Duley et al 2000). This would occur in regions of 
shocks and high UV excitation in the disk. The smaller PAHs 
may then be more easily 'levitated' above the plane. Jones et 
al. Q996) have also noted that PAHs and smaller particles can 
build up at the expense of larger ones in shocks or winds and 
that wind velocities as low as tens of km/s can modify the grain 
properties. The apparent dearth of /1850 yum emission in the 
halo (Sect. 14.2b supports this view. It would be interesting to 



obtain more sensitive A850 fim and Al. 2 mm data in order to 
compare the results with starburst galaxies, such as NGC 5775 
in which /i850 fim emission has been detected in the halo (Brar 
et al. 2003 1. Determining the admixture of different dust com- 
ponents as a function of height would provide important clues 
as to what mechanisms may be involved. Conditions in galaxy 
halos are quite different from those in the disk and may more 
closely resemble those seen in elliptical galaxies. 

Although there is a large optical ring around NGC 5907, as 
described in Sect. 11.11 there is no evidence that the high lati- 
tude PAH emission is directly related to it. The intersection of 
the optical ring with the disk of NGC 5907 occurs at two lo- 
cations, one of which (to the north) is not in the field of view 
of our maps and the other is near the southern tip of the major 
axis (coordinates given in Fig. II li . We see no evidence for en- 
hanced emission or extensions at the latter location. There does 
appear to be a disturbance in the Ha map of Rand et al. ( 1006) 
approximately at this location, however. 

6. Conclusions 

We have mapped the MIR emission of the low SFR galaxy, 
NGC 5907, in the 6.7W, 6.8N, 7.7N, 9.6N, 11. 3N, and 12W 
(LW2, LW5, LW6, LW7, LW8, and LW10, respectively) ISO 
bands, at high sensitivity and high resolution. Although we do 
not have spectra for NGC 5907, we have compared the line 
strengths and line ratios with those of other galaxies, especially 
M 83 in quiescent regions between the spiral arms. These com- 
parisons indicate that the excitation conditions are quite typi- 
cal of regions and/or galaxies with a low SFR. The contribu- 
tion from hot VSGs is likely no more than ~ 10% across these 
bands and is likely negligible in all bands except the broad, 
IRAS -equivalent 12W band. The MIR emission is dominated 
by PAH emission. 

In the disk, the high spatial resolution (1") broad band 12W 
emission is peaked towards the optical dust lane. The bands 
dominated by PAH emission (at 6" resolution) are observed 
all along the major axis of the galaxy and extend almost to 
the end of the optical disk (cf ~ 35 kpc in comparison to 41 
kpc) and possibly farther. The distribution of PAHs in the disk 
roughly follows that of the CO and not the broad HI distribu- 
tion. The PAH distribution shows an apparent correlation with 
the SCUBA /1850 fim distribution to within the positional er- 
rors of the data. The Al. 2 mm emission also follows that of the 
molecular gas. At the same time, there is some PAH emission 
(and also Al. 2 mm emission) at larger galactocentric radii than 
the CO, suggesting that PAH emission is more widespread than 
the molecular gas alone. It is likely that greater PAH emission is 
seen where there is more molecular gas because there are more 
UV photons available from SF in these regions and because 
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of the greater concentration of grains and PAHs in higher den- 
sity gas. A smaller fraction of the PAH emission will be com- 
ing from the more extended HI, excited by the general diffuse 
ISRF. 

An important finding of this work is the detection of PAH 
emission at high latitudes. The emission features show struc- 
ture in the narrow bands that isolate single PAH features, ap- 
pearing similar to features seen in Ha, radio continuum, and 
other bands in other edge-on galaxies (see, e.g. Lee et al. 2002 1. 
Individual features extend to latitudes of up to 6.5 kpc but 
a more typical scale height is 3.5 kpc. Previous high latitude 
emission on kpc scales has been seen only in the radio contin- 
uum component in this galaxy. Thus, cosmic rays and magnetic 
fields are known to exist in the halo, but other components have 
not (yet) been observed. Indeed, the PAHs appear to be selec- 
tively represented in the halo in comparison to large grains. 
We note also an increase in the Al 1.3//T7.7 PAH ratio with dis- 
tance from the major axis. Although NGC 5907 is a low SFR 
galaxy, there is some SF activity occurring, such as has been 
seen via some anomalous CO emission as well as a previously 
unreported off-nucleus radio continuum point source which is 
likely a supernova remnant or remnants. Thus, the halo PAHs 
may be transported to high latitudes with the help of coupling 
to magnetic fields (likely with gas as well), from regions of SF. 

Although NGC 5907 is surrounded by a faint stellar ring, 
there is no indication in these data that interaction with the ring 
plays a direct role in the presence of the high latitude emission. 
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